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Abstract 

We have investigated the invariant mass spectrum of dimuons collected by the CDF experiment 
during the 1992 — 1995 run of the Fermilab Tevatron collider to improve the limit on the existence 
of narrow resonances set by the experiments at the SPEAR e + e~ collider. In the mass range 
6.3 — 9.0 GeV/c 2 , we derive 90% upper credible limits to the ratio of the production cross section 
times muonic branching fraction of possible narrow resonances to that of the T(15) meson. In 
this mass range, the average limit varies from 1.7 to 0.5%. This limit is much worse at the mass of 
7.2 GeV/c 2 due to an excess of 250±61 events with a width consistent with the detector resolution. 
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I. INTRODUCTION 



Supersymmetric theories provide a mechanism to break the electroweak symmetry and 
to stabilize the large hierarchy between the Planck and the Fermi scales. Supersymmetry 
requires the existence of scalar partners to each standard model (SM) fermion, and spin— 1/2 
partners to the gauge and Higgs bosons. In particular, supersymmetry predicts the exis- 
tence of scalar quarks, i.e. particles that carry color, but no spin. Scalar quarks (squarks) 
have been searched for at current and past colliders, but charge— 1/3 squarks might have 
been overlooked for several reasons. Within the minimal supersymmetric standard model 
(MSSM), Ref. lj shows that the existence of a yet undetected charge— 1/3 scalar quark, 
lighter than the b quark, would require a lot of fine-tuning of the MSSM parameters, but at 
present cannot be ruled out by the electroweak precision data and the Higgs mass constraints 
from LEP2. Charge— 1/3 squarks would give a 2% contribution to the inclusive cross section 
for e + e~ — > hadrons; this contribution is comparable to the experimental error of the present 
measurements . Searches for narrow resonances at SPEAR have set upper bounds on 
Ti, the leptonic width of possible resonances, of 100 eV in the mass region 5.7 < E cm < 6.4 
GeV and of approximately 60 eV in the region 7.0 < E cm < 7.4 GeV [4]. In Ref. [51], the lep- 
tonic width of 1 bound states of such squarks has been evaluated using potential models 
of ordinary heavy quarks Q|. Because of the p-wave suppression of the fermion contribution 
to their decay width, the leptonic width is estimated to be approximately 18 (6) eV for a 
resonance with a 6 (10) GeV/c 2 mass. As noted in Ref. |J, for quarkonium masses above 
6 GeV/c 2 the width Ti is well below the experimental bounds, and scalar-quark resonances 
might have been missed. With this study, we investigate the region above 6.3 GeV/c 2 by 
using muon pairs with invariant mass between the ip and T mesons. The large statistics data 
set has been collected with the Collider Detector at Fermilab (CDF) during the 1992 — 1995 
Fermilab collider run. Section |H] describes the detector systems relevant to this analysis. 
Section ITTT1 reviews the expectation for bound states of charge —1/3 squarks. The data sam- 
ple is described in Sec. II V| and Section IVl presents additional selection criteria, tuned with 
T(1S) decays, which reduce the non-resonant background by a factor of three without losing 
more than 10% of the signal. In Sec. IVI[ we fit the dimuon invariant mass distribution and 
derive a 90% Bayesian upper limit on Ti as a function of the resonance mass. The shape of 
the invariant mass distribution is generally quite smooth and we improve the SPEAR limit 
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by an order of magnitude. An exception is the mass of 7.2 GeV/c 2 , at which the data can 
accommodate a Gaussian bump of 250 ± 61 events. In Sec. I VI II we explore the possibility of 
observing a real signal, we estimate its statistical significance, and we study its robustness 
by using a number of different kinematical selections. Our conclusions are summarized in 

Sec. EES] 

II. CDF DETECTOR 

CDF is a multipurpose detector, equipped with a charged particle spectrometer and a 
finely segmented calorimeter. In this section, we recall the detector components that are 
relevant to this analysis. The description of these subsystems can be found in Ref 
Two devices inside the 1.4 T solenoid are used for measuring the momentum of charged 
particles: the silicon vertex detector (SVX) and the central tracking chambers (CTC). The 
SVX consists of four concentric layers of silicon microstrip detectors surrounding the beam 
pipe. The CTC is a cylindrical drift chamber containing 84 sense wire layers grouped into 
nine alternating superlayers of axial and stereo wires. Electromagnetic (CEM) and hadronic 
(CHA) calorimeters surround the tracking volume and measure energy deposits over the 
pseudorapidity region \r)\ < 1. Muons are reconstructed by matching track segments in the 
drift chamber systems located outside the CHA (CMU, CMP, and CMX muon detectors, 
which cover the region \rj\ < 1) to the tracks of charged particles reconstructed in the 
CTC. The dimuon events used in this analysis were collected with a three-level trigger 
system. The first level required two charged tracks in the muon chambers. The second level 
trigger required that both muon tracks match a charged particle with transverse momentum 
Pt > 2.2 GeV/c as measured by a fast track processor (CFT). The third level software 
trigger requires that two charged CTC tracks, fully reconstructed in three dimensions, match 
track segments in the muon chambers and that the dimuon invariant mass is larger than 
2.8 GeV/c 2 . 

III. SEARCH FOR NARROW RESONANCES 

In this study, we search for narrow resonances 6, bound states of scalar quarks, in the 
dimuon invariant mass distribution between 6.3 and 9.0 GeV/c 2 . For a charge— 1/3 squark, 
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the muonic width, T^E — > of 2P resonances in this mass range has been evaluated 

in Ref. to be approximately 15 eV. In contrast, the T(IS') meson has a larger leptonic 
width of 1.32 keV. In analogy with the T(1S[) meson, the annihilation of a E state into 
hadrons is believed to proceed through gluons |5j, the dominant contribution coming from 
the minimum number of intermediate gluons (one-gluon is excluded by color conservation 
and two-gluons by C-conservation). In pp collisions quarkonia states are directly produced 
through subprocesses such as gg — ► g¥(lS) or gg — > g E. As the production subprocesses 
are directly related by crossing to the corresponding decay processes, the production cross 
section a is determined by the decay widths 8|, and approximately reads: 

a £ B(E - MM) - ( ~Ws) a r(is) S(T(15) - fifi) = R a T{ls) £?(T(15) - W ) (1) 

V m E ) 

where i? is the branching ratio of the muonic decay and = y-^^J prfrsy- F° r an 
£ particle in the mass region investigated by this study, R is approximately 2% when the 
width is evaluated using the standard potential of heavy-quark spectroscopy p . Using 
a different potential model and one-loop corrections to the static potential of the scalar 
quark- ant i- quark system, Reference [^| predicts leptonic widths that are a factor of three 
larger. In conclusion, R is expected to be between 2 and 6%. Since CDF has collected 
approximately 10 4 T(IS') mesons, the data could contain at least 200 events contributed 
by a hypothetical E meson on top of the smooth background due to Drell-Yan production, 
double semileptonic decays of cc and bb pairs, and fake muons produced by hadrons that 
mimic their signal. 



IV. DIMUON DATA SAMPLE 



The dimuon sample used in this analysis corresponds to approximately 110 pb 1 of data 
collected with the CDF detector during the 1992—1995 collider run. This data set has 
been usee, in Severn. CDF analy8 e 8 and is described in m ore detail in Re, Q. Tbe ennon 
identification is based on the three-dimensional matching of the track segment in the muon 
chambers with the track reconstructed in the CTC and on the energy deposited in the 



calorimeter tower in the muon path Q, 
Pt > 3 GeV/c using the same selection criteria of Ref. [1 
Pt < 3 GeV/c using the stricter cuts of the SLT algorithm 



111 ]. In this st udy we search for muons with 
we select muons with 2 < 
12 . 113I ] in order to reduce the 



misidentification background. We require that at least one of the muons is identified by 
both the CMP and CMU systems (j>t > 3 GeV/c and pseudorapidity \r]\ < 0.6). Additional 
muons are identified by either CMU or CMX system (p? > 2 GeV/c and \r]\ < 1.0). These 
selection criteria yield reconstructed muon pairs with rapidity \y\ < 1.0. We retain events 
which contain two and only two muons. The muon momentum is evaluated with a fit which 
constrains the track to originate from the beam line. The dimuon invariant mass is calculated 
using these momenta. This study uses opposite charge dimuon pairs; however, distributions 
for same charge dimuon pairs are also shown as a cross-check. 

The dimuon invariant mass distribution is shown in Fig. ^ The yield of ip mesons is 
much suppressed with respect to that of T mesons. Because of the muon selection criteria, 
at the J / ip mass the kinematic acceptance decreases rapidly with J / ip decreasing transverse 
momentum and vanishes at pr — 5 GeV/c. In contrast, for muon pairs with invariant masses 
larger than 6.3 GeV/c 2 , the kinematic acceptance does not depend on the dimuon transverse 
momentum. Therefore, we avoid the uncertainty of modeling an acceptance that depends on 
the production kinematics by limiting our study to the mass region above 6.3 GeV/c 2 . The 
dimuon invariant mass distribution at the T(IS') and in the region of interest for this study 
are shown in Fig. |21 The number of T(IS') mesons in the data is derived by fitting a first 
order polynomial and a Gaussian function to the invariant mass distribution in Fig. |2fa) 
with a binned maximum likelihood method. The best fit returns 9838 ± 141 Y(IS') mesons 
over a background of 5769 events in the region 9.3 < M At + At - < 9.55 GeV/c 2 . The fit also 
returns Mr(is) = 9439 ± 1 MeV/c 2 and a mass resolution au = 57 ± 1 MeV/c 2 . This mass 
resolution is well modeled by a simulation of the process pp — > T(1S)X. The simulation 
event generator produces T(IS') mesons with the transverse momentum distribution of the 
data Q and a flat rapidity distribution for | y | < 1. The generated events are proeessed with 
the CDF detector simulation 1 (qfl) described in detail in Ref. [15j. Events are then required 
to pass the same selection and reconstruction criteria imposed on the data. This simulation 
predicts a mass resolution of 40 MeV/c 2 for E states with a mass around 7.5 GeV/c 2 . 



The simulation includes correction factors for the efficiency of the three-level trigger system, the efficiency 
for reconstructing tracks in the CTC and the different muon systems, effects due to instantaneous lumi- 
nosity and to internal radiation from muons. These correction factors are parametrizations based on the 
data QE1- 
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FIG. 1: Invariant mass distribution of all muon pairs. OS indicate dimuons with opposite sign 
charge used in this analysis. Dimuons with same sign charge (SS) are also shown. 
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FIG. 2: Invariant mass distributions of muon pairs at the T(1S) (a) and in the region of interest 
for this study (b). The solid line represents the fit used to estimate the number of T(IS') mesons. 

V. BACKGROUND REDUCTION 

As outlined in Sec. IIII| the hypothetical signal of an 6 resonance of mass ~ 7.5 GeV/c 2 
is expected to be at least 2% of the T(IS') yield, i.e. about 200 events. Given the detector 
invariant mass resolution, the 200 events have to be integrated in a region of 150 MeV/c 2 
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(±2<tm)- As shown in Fig. Efb), three 50 MeV/c 2 bins centered around this mass contain 
approximately 8000 events. This background can be largely suppressed because it is mostly 
contributed by bb and cc production. The measurement of the time-integrated B° — B° 
mixing probability, reported in Ref . , also makes use of this data sample. From that study 
we estimate that approximately 75% of the muon pairs arise from heavy flavor production. 
We use two intuitive criteria to reject dimuons arising from the decay of hadrons with heavy 
flavor: 

1. An isolation requirement. The isolation, /, is defined as the scalar sum of the transverse 
momenta of all the tracks in a cone of radius R = y/5<p 2 + 5r]' 2 = 0.4 around the muon 
direction. We require that both muons have isolation / < 4 GeV/c. 

2. A promptness requirement. In contrast to b and c-hadrons, the T(IS') and 6 mesons 
have negligible lifetime. We select prompt muons by requiring the sum of the impact 
parameter significance of both muons, Si P , to be less than 4 2 . The impact parameter 
significance is estimated for muons with tracks reconstructed in the microvertex silicon 
detector (SVX), otherwise is set to zero in order not to lose events. 

The values of these cuts have been determined at the T(IS') mass (see Fig. 0J). As shown 
in Table HJ these cuts reduce the background by more than a factor of three while retaining 
more than 90% of the T(15') signal. 



2 The track impact parameter d is the distance of closest approach to the event primary vertex in the plane 
transverse to the beam line. The significance is defined as d/<Jd- The event primary vertex is determined 
as in the study in Ref. Q. 
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FIG. 3: Invariant mass distributions of opposite (•) and same sign charge (dot-dashed) muon pairs 
in the T (IS) mass range after the isolation (a) and impact parameter cuts (b). The solid line is 
a fit using a Gaussian plus a first order polynomial functions. 



TABLE I: Numbers of T(15) mesons and underlying-background events for different analysis cuts. 
Rates are evaluated by fitting the data with a Gaussian plus a first order polynomial function. The 
background is integrated over the mass region between 9.3 and 9.55 GeV/c 2 . The last cut, QC, is 
used in the 6 search. Efficiencies are calculated with respect to the number of T(15) candidates 

in the first row. 

Cut T(15) candidates Background Efficiency(%) 

None 9838 ± 141 5769 

I 9821 ± 129 3345 99.8 

QC = I + s ip 9064 ± 118 1842 92.1 
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VI. ESTIMATE OF THE UPPER LIMIT ON if 



Figure 0] shows the invariant mass distribution of muon pairs in the region between 6.3 
and 9.0 GeV/c 2 after the isolation and impact parameter significance cuts. 

We use a binned maximum likelihood method to fit the mass spectrum in Fig. 0] with a 
fourth order polynomial, which serves the purpose of modeling a smooth background, plus 
a Gaussian function, which searches for narrow resonances. We perform 54 fits, in which we 
constrain the Gaussian peak to the center of each of the 54 mass bins of Fig. HJ in each fit, 
we force the Gaussian width to the simulated resolution of the detector for that mass. For 
each mass bin, we use the integral of the Gaussian function and its error returned by the 
best fit to derive N u i, the 90% credibility upper limit to the number of events contributed by 
a narrow resonance centered in that mass bin 3 . We evaluate the ratio of the geometric and 
kinematic acceptance for an £ resonance to that for the T(1S) meson with the simulation 
described at the end of Sec. IIVI In the event generator, an unpolarized resonance is produced 
with a flat rapidity distribution (\y\ < 1) and a transverse momentum distribution that is 
rescaled from that of the T(IS') data so that < > / < Pt > = m £/ m r(is)- For J/ifj 
mesons, this rescaling procedure predicts a da/dpr distribution that decreases more rapidly 
with increasing momenta than the distribution of the data. However, a poor modeling of the 
transverse momentum distribution is not a cause of error because the kinematic acceptance 
does not depend on the £ transverse momentum. As a cross-check, we also generated £ 
resonances using the shape of the transverse momentum distribution of the T(IS') data 
and verified that they return consistent acceptance values. The geometric and kinematic 
acceptance increases from 8.3% at 6.3 GeV/c 2 to 11.2% at 9.0 GeV/c 2 (it is 11.7% at T{IS) 
mass). The parametrized correction factors for the trigger and reconstruction efficiencies, 
discussed at the end of Sec. IIVI depend little on the muon px- Including these effects in 
the simulation, the acceptance is 4.7% at 6.3 GeV/c 2 and 6.3% at 9.0 GeV/c 2 (it is 6.6% 
at T(15') mass). The parametrized correction factors to the acceptance have approximately 
a 6% uncertainty. However, in the mass range considered in this study, the ratio of the 
simulated £ to T(IS') acceptances is not affected by this uncertainty. 

3 The integral of the fit likelihood from N u i to infinity is 10% of the integral of the fit likelihood from to 
infinity. The justification for this procedure is Bayesian with a prior that is zero for negative resonance 
cross sections and flat for positive ones. 
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The ratio of N u i to the number of observed T(1S) mesons, corrected for the relative 
acceptance, provides the 90% credibility upper limit to erg- B(£ — ► fifi) / arns) B(T(1S) — > 
fi/j,) shown in Fig. El Figure El shows the 90% credibility upper limit to derived using 
equation (1). We note that equation (1) tends to underestimate the £ production cross 
section. One can verify this by predicting the J /if) production cross section from that of 
the T(IS') meson that is about two order of magnitude smaller. The T(IS') integrated cross 
section for \y\ < 0.6 is measured to be 34.6 ± 2.5 nb 4 . The analogous cross section for 
the J /if) meson is 3060 ± 460 nb 5 , whereas equation (1) predicts 1630 ± 116 nb. If the £ 
production cross sections predicted by eqn. (1) were also a factor of two smaller than the 
data, the vf limits set by our study would be a factor of two smaller than those indicated 
in Fig. El Figure [7| shows the distribution of iVg-, the number of events attributed to a 
narrow resonance, divided by the error returned by the best fit for the 54 considered 
mass bins. With the exception of a point at 4.1 a, this distribution is consistent with a 
Gaussian function of unit width. Therefore, it seems fair to assume that the distribution of 
the 90% upper limits in Fig. El is statistically consistent with the average upper limit = 8 
eV that corresponds to those cases with Ng/<jN = in Fig. The 4.1 a fluctuation occurs 
at the mass of 7.25 GeV/c 2 . 



4 Average of the measurements in Refa, |14| . 

5 We use the result reported in Ref. |l7j ]. The measurement is performed at \fs = 1.96 TeV, and we rescale 
the published value by the 10% expected increase of the cross section from 1.8 to 1.96 TeV 0|- We also 
extrapolate the result to p™ m = GeV/c assuming that the fraction of prompt J/ip remains constant for 
transverse momenta smaller than 1.5 GeV/c. 
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FIG. 4: Invariant mass distribution of opposite (•) and same sign charge (dot-dashed) muon pairs 
which pass the isolation and impact parameter cuts. The solid line represents the fit described in 
Sec. EIU 
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VII. STUDY OF THE 7.25 GeV/c 2 FLUCTUATION 

In this section we explore the possibility that the anomalously large upper limit at 
7.25 GeV/c 2 is due to a real signal. A fit, which uses a Gaussian function with a fixed 
38 MeV/c 2 resolution, returns M £ = 7.22 ± 0.01 GeV/c 2 and a signal of 250 ± 61 events 
over a smooth background of 3355 events, extrapolated in the 150 MeV/c 2 region between 
7.15 and 7.3 GeV/c 2 (see Fig. HJ). The probability that 3355 background events fluctuate 
to no less than 3605 is 8 x 10~ 6 (4.3 a). Since the mass range examined in Fig. 0] includes 
52 almost independent combinations of three consecutive 50 MeV/c 2 bins, the probabil- 
ity of obtaining an equal or larger statistical fluctuation in the inspected mass window is 
approximately 4.1 x 10~ 4 (3.5 a). 

According to the simulation the acceptance for a 7.2 GeV/c 2 resonance relative to that 
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for the T(IS') meson is = 0.78 x Arns)- It follows that 

a £ B(£ -> fifi) = (3.6 ± 0.9) x 10~ 2 x cr r{is) B(T(1S) -> ////) 

We note that this value is in agreement with the theoretical expectation for a bound state 
of charge— 1/3 squarks. 

We have investigated three additional selection cuts that reduce the number of events 
by more than a factor of two and compare the effect of these cuts on the number of £ and 
T(15 l ) candidates: 

1. | cos(#*)| < 0.4, where 9* is the polar angle between the /i + and £ directions in the £ 
center-of-mass system. This cut reduces the number of events by simply selecting a 
particular sector of the phase space. 

2. y^p r < 40 GeV/c and < 30, where ^2p l T is the scalar sum of the transverse 

i j i 

momentum of all tracks i originating from the same vertex as the muon pair and y~^ p 

3 

is the sum of the impact parameter significance of all tracks j not used to define the 
primary vertex of the event [18j. This cut is intended to further suppress the heavy 
flavor background by rejecting events in which the muon pair carries a small fraction of 
the total transverse momentum or is produced in association with additional long-lived 
particles. 

3. cut #1 + cut #2 

The effect of these cuts is shown in Figs. |S1 and El and is compared to the result for the 
T(IS') meson in Table ITT1 
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TABLE II: Yield of £ and Y(1«S) candidates for different analysis cuts. The underlying backgrounds 
are fitted with polynomial functions and integrated over the mass ranges 7.15 — 7.3 and 9.3 — 
9.55 GeV/c 2 . 



Cut £ Background T(1S) Background 

QC 249.7 ±60.9 3355.0 9064 ± 118 1842 

#1 160.5 ±41.8 1508.0 3910 ± 90 611 

# 2 206.2 ± 57.0 2948.0 8667 ± 136 1587 

#3 144.5 ±39.0 1311.0 3699 ± 87 561 
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FIG. 9: Invariant mass distributions of opposite sign charge muon pairs after applying the third 
cut described in Table ITT1 
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VIII. CONCLUSIONS 



We have investigated the invariant mass spectrum of dimuons collected by the CDF 
experiment at the Tevatron collider to improve the limit to the existence of narrow 
resonances set by the experiments at the SPEAR e + e~ collider. In the mass range 
6.3 — 9.0 GeV/c 2 , we derive 90% upper credible limits to the ratio of the production cross 
section times muonic branching fraction of possible narrow resonances to that of the T(IS') 
meson. In this mass range, the average limit varies from 1.7 to 0.5%. Assuming that 
ag = f7x(i5) x (m T (!5)/m^) 3 x /Tj^ is \ these limits correspond to an average 90% upper 
credible limit of 8 eV to the leptonic width of possible resonances. An exception is the 
mass region around 7.2 GeV/c 2 where we observe a bump of 250 ± 61 events with a width 
consistent with the detector resolution. The size of the excess is consistent with the theoret- 
ical expectation for the production of a 1 p-wave resonance but its statistical significance 
(3.5 a) is not sufficient to claim the discovery of a new particle. 
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